Effective triplet interactions in nematic colloids 
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Three-body effective interactions emerging between parallel cylindrical rods immersed in a ne- 
matic liquid crystals are calculated within the Landau-de Gennes free energy description. Collinear, 
equilateral and midplane configurations of the three colloidal particles are considered. In the last 
two cases the effective triplet interaction is of the same magnitude and range as the pair one. 
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I. INTRODUCTION 

In colloidal suspensions knowledge of particle interac- 
tions is the key to understanding their thermodynamic 
properties. Typically it is assumed that the main contri- 
bution comes from the pairwise interactions between the 
particles. This simplifies the corresponding equations of 
motion which can then be solved using efficient numeri- 
cal algorithms, such as molecular dynamics. However, in 
many systems the governing physical equations are non- 
linear and this, in principle, can result in a modification 
of the interaction between two particles by the presence 
of other ones, meaning that the total interaction energy 
is no longer pairwise additive, and additional many-body 
potentials appear. The importance of many-body inter- 
actions has been recognized in a number of systems, from 
electron screening in metals jl| , catalysis and island for- 
mation on surfaces [1] , to charged colloids [p, 0, H, S 0] j 
binary mixtures 0] , colloidal crystals 0, [lOI | and simple 
liquids [IIGJ]. 

Attempts to characterize the solvent-mediated inter- 
actions have been made in case of nematic colloids [H, 
EH EH > where colloidal particles experience long-ranged 
interactions mediated by the distortions of the nematic 
director around the particles [lfl E3, EH- These effec- 
tive interactions result in particle clustering and self- 
organization Due to the presence of topological 
defects [2(], |21|. l22t [231. l2~j ] the interparticle interaction 
is inherently anisotropic: depending on the symmetry of 
the director field around the particles either dipole-dipole 
or quadrupolar-type interactions can be observed [l7ll2{j ] . 

The long-ranged nature of interactions mediated by 
the nematic director can result in a substantial many- 
body contribution to the interaction energy already for 
small concentrations of colloids. Indeed, recent theoret- 
ical studies of Aimers [26| and trimers [23] dispersed in 
a nematic host have shown that the three-body interac- 
tions can be of the same importance as two-body ones. 

In spite of this fact, most research on nematic colloids 
still relies on pairwise interaction potentials between the 
particles [H, H3, EL H3| ■ Moreover, often the direc- 
tor distribution around N particles is approximated by 



a superposition of one-particle solutions of a correspond- 
ing linearized theory (so-called Nicolson approximation), 
which yields vanishing many-body potentials. In addi- 
tion, a linear superposition of one-particle solutions does 
not satisfy boundary conditions at all particle surfaces. 

In this paper we would like to address the role of 
three-body interactions for long cylinders immersed in 
a nematic mesophase. This system can serve as a two- 
dimensional representation for nanoparticles in a nematic 
host. It is also a realistic model for carbon nanotubes- 
nematic liquid crystal composites (discotic carbonaceous 
mesophases), which have recently been studied theoreti- 
cally [26[ and experimentally [33| . 



II. MANY-BODY POTENTIALS AND 
NEMATIC FREE ENERGY 

Due to the anisotropic nature of the solvent, an analyt- 
ical treatment of particle interactions is practically im- 
possible already at a pairwise level. We therefore choose 
here a numerical approach: equations for the nematic 
order tensor are solved by minimizing the free energy 
numerically using finite elements with adaptive meshes 
|34| . For a system of N colloids at positions we in- 
troduce an excess free energy Fn (defined below) over 
the free energy of a uniaxial nematic. Fn is then decom- 
posed into effective n-body potentials F^ n \ with n < N, 
as following Q 

JV N 

F N = NF^+J2F (2) (hj)+ E F (3) (i>3,k) + ..., (1) 

i<j i<j<k 

where we introduce the short-hand notation for the cen- 
ter of mass positions (rj, ...) = (£,...). The n-body po- 
tential is defined in the system with n colloidal particles. 
i^ 1 ) = F\ is the excess free-energy coming from one iso- 
lated colloid. The effective pair interaction between two 
colloids follows from Eq. (TTJ) when N = 2 

FW(l,2)=F 2 {l,2)-2F 1 . (2) 
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The triplet potential (1, 2, 3) is defined in the system 
of three colloids, N = 3 

(i, 2 , 3) = F 3 (l, 2, 3) - ]T F< 2 > - 3F X . (3) 

By construction F^(l, 2, 3) tends to zero whenever one 
of the Tij = \r,i — rj\ — > oo. ^^(1,2,3) is the quantity 
we are interested here. 

The nematic host is modelled by the Landau-de 
Gennes free energy density [35j 

/ = a(T — T*)QijQji — bQijQjkQki + c(QijQji) 2 

~T" ~^~Qij,kQij.k ~t~ .jQik,k j (4) 

where is the order parameter tensor, summation over 
repeated indices is implied and the comma indicates the 
spatial derivative. The positive constants a, 6, c are as- 
sumed to be temperature independent, and T* is the 
supercooling temperature of the isotropic phase. The 
constants L\ and L2 are related to the Frank-Oseen 
elastic constants Ku = K 33 — 9<2f(Zi + L 2 /2)/2, 
K22 — 9Q 2 Li/2, and Qb is the bulk nematic order pa- 
rameter. We introduce the dimensionless temperature 
r = 2Aca(T - T*)/b 2 . The bulk nematic phase is stable 
for r < tni = 1 with a degree of orientational order given 
by Q b = 6(1 + y/1 - 8r/9)/8c; Q b {r > 1) = 0. 

For the constants entering the free energy density 
@ we use typical values for a nematic compound 5CB 
[SI [37[: a = 0.044 x 10 6 J/m 3 K, b = 0.816 x 10 6 J/m 3 , 
c = 0.45 x 10 6 J/m 3 , Li = 6 x 10~ 12 J/m, L 2 = 12 x 
10~ 12 J/m, T* = 307 K. The nematic-isotropic transi- 
tion temperature for 5CB is Tni = 308.5 K. We choose 
T = 308.3 K, corresponding to Qb rj 0.31, close to the NI 
transition point, at which flocculation of colloidal parti- 
cles, driven by attractive interactions, normally occurs. 

1/2 

The nematic coherence length is £ = (24Lic/6 2 ) w 
10 nm; it sets the smallest length-scale in our system. 

Three identical colloidal particles, each of which we 
take to be a long cylinder of radius R with the symmetry 
axis parallel to the z axis, are immersed into a nematic 
liquid crystal. The order parameter at the system bound- 
aries is fixed to the bulk order parameter of the nematic 
phase; the director orientation at the box boundaries is 
fixed along the y axis. In this geometry the colloidal par- 
ticles are orthogonal to the director field far from them 
(the other geometry is when the rods are aligned with 
the nematic) . The problem of equilibrium orientation is 
discussed in Ref. [38j |. 

The order parameter at the particle surfaces is fixed to 
the bulk order parameter. We consider fixed homcotropic 
boundary conditions at the particle surfaces, i. e. the 
director orientation is fixed along the surface normals. 

Assuming translational symmetry of the system along 
the z axis, the equilibrium distribution of the order pa- 
rameter tensor Qij minimizes the free energy functional 
F = J fdxdy, where / is the free energy density ([!]). 




FIG. 1: (Color on-line) Order parameter and director dis- 
tribution around colloidal particles with R/£ = 10. Left: 
»*23 = 3i?, ri2 = 3.5R; right: rn = 3i?, ris = 2.5R. Di- 
rector orientation far from the particles is parallel to the line 
connecting particle centers. Note that only a small part next 
to the particles is shown, the actual box size is 40i? x 40i?. 

The functional F is minimized numerically, using finite 
elements with adaptive meshing. The two-dimensional 
domain of integration of the free energy density / is tri- 
angulated and the functions Qij are linearly interpolated 
within each triangle. Then the excess free energy Fn is 
calculated as 

F N = L z J(f- f b )dxdy, (5) 

where /{, = f{Qb) and L z is the length of colloidal par- 
ticles. Here we will only talk about effective interactions 
per length, F N = F N /L Z , F^l = F (n 1/L z . It will be 
understood that such scaling has been carried out, and 
we shall omit the overbars in the text below. 



III. RESULTS 

A. Collinear geometry 

In this case the three colloids are positioned on a line 
parallel to the nematic director n, along the y axis. Typ- 
ical director orientation and order parameter maps are 
shown in Fig. ([1]). Each particle is accompanied by two 
defect lines of strength —1/2, and the orientation of each 
defect pair is tilted with respect to the (fixed) director 
orientation far from the particles, Fig. (QJ,). Note that 
if a single colloidal particle is immersed in the nematic 
host, the accompanying defects will be positioned on a 
line perpendicular to the director orientation. The tilt 
angle (which breaks the mirror symmetry of the system) 
is a pure two-body effect and is already present in a sys- 
tem of two colloids [39]. At small separations, the three 
particles in a collinear geometry can serve as a model 
for a rod for which symmetry breaking has also been ob- 
served [40]. When the distance between the particles (1) 
and (2) is decreased (see Fig. (JTJd)), the orientation of the 
defects next to these particles tilts even more. However, 
this does not affect the order parameter and the director 
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FIG. 2: (Color on-line) Collinear geometry. Three colloids 
are positioned on the line which is parallel to the director 
n. Particle radius 7?/£ = 10. (a) Excess free energy F3 as 
a function of the separation between the particles (1) and 
(2). The position of the particle (3) is fixed with respect to 
the particle (2) at 37?, 4R, and 57?. (b) Pair potentials for 
r23 = SR. (c) Three-body contribution to F3; solid, dashed 
and dash-dotted lines correspond to 7-23 = 37?, 47? and 57?, 
respectively. Inset shows the ratio between the three- and 
two-body contributions and a power law fit of the decaying 
tail for r23 = SR. 



distribution around the third particle, i. e. one might 
anticipate that in a collinear geometry the contribution 
of the three-body potential to the total interaction en- 
ergy is rather small, in line with the results of Guzman 
et. al. H3|. 

To asses this contribution, the total interaction energy 
F3 as well as the pair and the triplet contribu- 
tions to it are shown in Fig. @ . Here we fix the distance 
between the particles (2) and (3) at 3i?, 4i?, and 57? 
and study the interaction of these two particles with the 
particle (1). The dependence of F$ on separation 7*12 is 
similar for all three curves: particle (1) is repelled from 
the particles (2) and (3) at large distances and there is a 
barrier which does not allow particle (1) to form a sta- 
ble aggregate. The height of this barrier decreases with 
the decrease of the distance between the particles (2) and 
(3). Figure (fJb) shows that the major contribution to the 
total interaction energy comes from the pairwise interac- 
tion of particles (1) and (2), especially at large separa- 
tions of particles (2) and (3). However, for r^j, — 37? we 
have significant contribution of the three-body potential, 
which is shown in Fig. (T3J:). In fact, this contribution is 
responsible for the decrease of the repulsive barrier and 
will lead to more efficient clustering of the colloids. For 
r 23 iS 37? the pairs of defects around particles (2) and (3) 
start to tilt. The tilt angle increases with the decrease 
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FIG. 3: (Color on-line) Collinear geometry, dependence on 
the tilt angle with respect to the nematic director: (a) ex- 
cess free energy 7<3; (b) two-body potentials; (c) three-body 
potential, solid, dashed and dash-dotted lines correspond to 
ri2 = 4.427?, 5.667? and 6.907?, respectively. Particle radius 
R/Z = 10. 



of , affecting the effective interaction between the pair 
and the third colloid [4lj | . 

One can also see that in this geometry the three-body 
potentials always provide an attractive contribution to 
the total interaction energy. To emphasize the role of 
the three-body term, its ratio to the sum of two-body 
potentials is shown in the inset of Fig. (J2fc) for = 37?. 
It clearly shows that at small separations, when rear- 
rangement of the defects plays a major role, the three- 
body term is dominating and therefore, cannot be ne- 
glected. At larger distances the role of the three-body 
term becomes less important decaying as a power law 

We also study how the two contributions change under 
the director rotation. Figure ([3]) shows the tilt angle de- 
pendence of the excess free energy F3 , as well as two- and 
three-body contributions to it. It is interesting that our 
calculations predict an equilibrium configuration with a 
tilt angle about 70° with respect to the nematic director, 
i. e. there is also a torque acting on the line of colloids. 
The tilt is due to the pairwise interactions: the three- 
body contribution, in fact, has a maximum at about 60° 
and favors almost collinear with the director alignment 
(minimum at about 10°). 



B. Equilateral triangle geometry 

The order parameter and the effective interaction po- 
tentials for the equilateral triangle geometry are pre- 
sented in Fig. (H|) and Fig. ©, respectively. Already 
the director and order parameter profiles, depicted in 
Fig. (f?|) , show that the three-body potential is significant 
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FIG. 4: (Color on-line) Order parameter and director dis- 
tribution for the case of equilateral triangle geometry, ri2 = 
ris, — T2z = 2.8R, R/£ = 10. (a) three particles; (b) the same 
geometry, but the particle (3) is removed; (c) the particle (1) 
is removed. The director far from the particles is parallel to 
the line connecting the centers of particles (1) and (2). Only 
a small area next to the particles is shown. 




FIG. 5: (Color on-line) Equilateral triangle geometry, R/£ — 
10. (a) total interaction energy F3; (b) two-body contribu- 
tions; (c) three-body contribution. Inset shows the ratio be- 
tween the three-and two-body contributions and a power law 
fit (r~ 3 ' 5 ) of the tail. 



in this case: the position of the defects and the director 
orientation around three colloids (see Fig. cannot 
be obtained as a superposition of the corresponding two- 
particles profiles (see Fig. ((4t>) and Fig. ((It)), even at 
large separations. Since in the equilateral triangle ge- 
ometry defect structures are less separated than in the 
collinear one, a much stronger three-body contribution is 
expected in this case. 

In Fig. ((St) F3 as a function of the triangle size r is 
shown, demonstrating attraction of the particles at large 
r. Pairwise contributions show that the particles (1) and 
(2) repel each other at large distances, while pairs (2)-(3) 
and (l)-(3) attract. The three-body contribution is al- 
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FIG. 6: (Color on-line) Equilateral triangle geometry, tilt an- 
gle dependence of the effective interaction energies. R/£ = 10. 
Bottom panel: solid, dashed and dash-dotted lines correspond 
to r = 4A2R, 5.66R and 6.90R, respectively. 



most always positive, and has a minimum around r = 3i?. 
Contrary to the collinear geometry, the ratio F^ 3 ' / F^ 2 > , 
which is shown in the inset of Fig. ([5]), does not vanish at 
large r. Moreover, it converges to a constant (negative) 
value. This means that at large distances the three-body 
potential decays with the same power law as the two- 
body ones. Therefore, it cannot be neglected even at 
large separations between the colloidal particles. 

The effective interaction potentials reveal rather weak 
dependence on the tilt angle, which is shown in Fig. [6] 
The peaks around 30° are related to the defect rearrange- 
ment around the particles. 



C. Midplane geometry 

In this geometry the position of the particle (3) is fixed 
at a distance X3 from the line connecting particles (1) and 
(2) (the geometry is depicted in Fig. 0). The addition 
of the third particle changes the repulsive interaction be- 
tween particles (1) and (2) to the attractive one, favoring 
their flocculation. The F^ s '/F^ 2 ' ratio shows the same 
tendency as in the equilateral triangle case: it does not 
reach zero at large r 12 , i. e. the three-body potential 
is long-ranged and provides about 10% correction to the 
sum of pair potentials at all distances. 



D. Four-body interactions 

To get the feeling of the next (four-body) term in the 
expansion ([I]) we considered a system of four colloidal 
particles positioned at the vertices of a square. Two-, 
three-, and four-body potentials are shown in Fig. ©. 
The inset shows the ratio of the total three- and four- 
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FIG. 7: (Color on-line) Midplane geometry, R/£ = 10: (a) 
total interaction energy; (b) pairwise contributions; (c) three- 
body contribution; solid, dashed and dash-dotted lines corre- 
spond to xz = 2.5R, 3.5R and 5R, respectively. Inset shows 
the ratio between the three-and two-body contributions for 
x-i = 2.5R. 
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FIG. 8: (Color on-line) Effective interaction potentials for 
four colloidal particles in a square geometry, with R/£ — 10. 
(a) pairwise contributions; (b) three-body contributions; (c) 
four-body contribution. Inset shows the ratio between the 
three-, four- and two-body contributions with the power-law 
fits. 



F 



(2) 



F( 2 )(l,2) + F( 2 )(2,3) + F( 2 )(l,3) 



F (3) = 4F (3) (1,2,3). 

Fitting the ratios / F^ (i = 3, 4) to the power law 
(ai/x) li + bi in the region r G [6R..10R] yields 73 = 5.5, 
63 = 0.18 and 74 = 3.6, 64 = —0.01. Again, at short dis- 
tances the contribution of the three-body potentials to 
F 4 is comparable (and at small distances even bigger) to 
the contribution of the two-body potentials. Moreover, 
the fit shows that both contributions have the same range 
and that three-body potentials amount up to 10% of the 
two-body terms even at large separations. The four-body 
term is slightly smaller than the three-body one. How- 
ever, the slow convergence of the expansion ([1]) does not 
allow to make a quantitative estimate of a many-body 
contribution that can be neglected. 



IV. CONCLUSIONS 

Our calculations indicate that a collinear array of col- 
loids (in our case modeled as infinitely long rods) is not 
as stable as the equilateral triangle configuration. The 
most stable equilateral triangle configuration is when one 
of the sides of the triangle is parallel to the nematic di- 
rector. 

Russ at. el. Q have shown that the triplet interac- 
tion is a function of the summed distance r\2 + r23 + r^x 
between the three particles only. However, this empiri- 
cal observation does not hold in our case. The possible 
reasons are: anisotropic nature of the interactions and 
presence of topological defects, both of which are absent 
in a system of charged colloids. 

The main message of this work is that in nematic 
colloids many-body effects cannot be neglected. At 
short distances rearrangements of accompanying defects 
strongly contribute to the three- and four- body terms in 
the free energy expansion. In some geometries (e. g. equi- 
lateral triangle) the exponents of the power law decay of 
the two- and three- body contributions with interparticle 
distance are the same, i. e. the higher order contributions 
cannot be neglected even at large separations. 
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